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Abstract. A recent study by Yoon & Langer 1 2004a l indicated that the inner cores of rapidly accreting (M > 10 7 M G yr 1 ) CO 
[ — ' white dwarfs may rotate differentially, with a shear rate near the threshold value for the onset of the dynamical shear instability. 

Such differentially rotating white dwarfs obtain critical masses for thermonuclear explosion or electron-capture induced collapse 
which significantly exceed the canonical Chandrasekhar limit. Here, we construct two-dimensional differentially rotating white 
dwarf models with rotation laws resembling those of the one-dimensional models of Yoon & Langer 1 2004a l. We derive analytic 
relations between the white dwarf mass, its angular momentum, and its rotational-, gravitational- and binding energy. We show 
that these relations are applicable for a wide range of angular velocity profiles, including solid body rotation. Taken at a 
central density of 2 x 10 9 g cirT 3 they specify initial models for the thermonuclear explosion of rotating CO white dwarfs. At 
p c = 10'° g cirr 3 and 4 x 10 9 g cm" 3 , they give criteria for the electron-capture induced collapse of rotating CO and ONeMg 
' white dwarfs, respectively. We demonstrate that pre-explosion and pre-collapse conditions of both, rigidly and differentially 

rotating white dwarfs are well established by the present work, which may facilitate future multi-dimensional simulations of 
Type la supernova explosions and studies of the formation of millisecond pulsars and gamma-ray bursts from collapsing white 
dwarfs. Our results lead us to suggest various possible evolutionary scenarios for progenitors of Type la supernovae, leading to 
a new paradigm of a variable mass of exploding white dwarfs, at values well above the classical Chandrasekhar mass. Based on 
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our 2D-models, we argue for the supernova peak brightness being proportional to the white dwarf mass, which could explain 
various aspects of the diversity of Type la supernovae, such as their variation in brightness, the dependence of their mean 
luminosity on the host galaxy type, and the weak correlation between ejecta velocity and peak brightness. 
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' & 3 genitor scenarios, such as the sub-Chandrasekhar Mass sce- 

spun up to critical rotation, which influences significantly their . ^ , , . . , , , , 
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influence both, the SNe la progenitor evolution and the ther- 
monuclear explosion of the white dwarf. 

During the last two years, various systematic investiga- 
tions of the effects of the white dwarf spin-up due to accre- 
tion in the context of SNe la progenitor evolution have emerged 
(Tiersanti l2003al I2003bl Uenishi l2U0H Yoon & Langer|l004a 
I2004M Yoon et al. 120041 Saio & Nomoto l2004l . These results 
may demand surprising modifications of the standard picture 
of SNe la progenitors. 

Yoon, Langer & Scheithauer (2004i showed that the shell 
burning in rapidly accreting white dwarfs is stabilized by 
the accretion induced spin-up, as the shell source is widened 
through rotationally induced chemical mixing and decom- 
pressed by the centrifugal support (cf. Yoon, Langer & van der 
Sluys 2004 1. This result supports the key assumption of the 
Single Degenerate scenario that white dwarfs can efficiently 
grow in mass by shell burning, which has been challenged 
by the finding of generally unstable helium shell burning and 
subsequent mass loss in non-rotating white dwarfs (Cassisi et 
al. [TW8l Kato & Hachisu l2004l 

In the sub-Chandrasekhar Mass scenario, the detonation 
of a thick helium layer triggers the explosion of a sub- 
Chandrasekhar mass CO core (e.g. Woosley & Weaver 1994 1. 
Here, Yoon & Langer ( 2004b I showed that as the white dwarf 
is spun up by accretion at the moderate rates required to build 
up the degenerate helium layer, the latter is efficiently heated 
by the dissipation of rotational energy. Consequently, helium 
ignites much earlier and under less degenerate conditions than 
in the corresponding non-rotating case, which avoids a super- 
nova explosion altogether. 

Most intriguingly, a study by Yoon & Langer (2004a, here- 
after Paper I) indicated that accreting white dwarfs with M > 
10 -7 M Q yr~' may rotate differentially throughout. They found 
that, although the strongly degenerate inner core is suscepti- 
ble to the dynamical shear instability (DSI), other instabilities 
such as Eddington-Sweet circulations, the secular shear insta- 
bility and the Goldreich, Schubert and Fricke instability are 
relevant for angular momentum redistribution at the required 
time scales (~ 10 6 yr) only in the non-degenerate outer layers. 
Accordingly, the shear rate in the inner core of the white dwarf 
corresponds to the shear rate near the threshold value for the 
onset of the DSI (Paper I). Since differentially rotating white 
dwarf can be dynamically stable even up to 4.0 M , the masses 
of exploding white dwarfs could be significantly larger than 
the canonical Chandrasekhar limit of ~1.4 M (cf. Paper I).. 
This may require to supersede the long-standing paradigm of 
Chandrasekhar-mass white dwarfs as SN la progenitors with 
a new paradigm of SN la progenitors as white dwarfs with a 
variety of masses at super-Chandrasekhar mass values. 

The possibility of super-Chandrasekhar mass white dwarfs 
is not new; already in the 1960ies, e.g. Ostriker & 
Bodenheimer (1968 1 constructed differentially rotating 2-D 
white dwarf models of up to 4.1M . However, the internal an- 
gular momentum distribution in these models had to be as- 
sumed and was chosen such that the spin rate decreases out- 
ward — conversely to the result of evolutionary considera- 
tions (Paper I). More recently, Uenishi et al. (|2003l provided 
2-D white dwarf models, assuming a rotation law which in- 



troduces a discontinuity between the non-rotating inner core 
and the rigidly rotating spun-up envelope. However, due to this 
assumption, the maximum possible white dwarf mass became 
1 .47M , which strongly restricts the diversity of possible ex- 
plosion masses. In the present paper, we construct 2-D white 
dwarf models with rotation laws that mimic the 1-D results of 
Paper I and thus follow from calculating the evolution of the 
internal white dwarf angular momentum profile throughout the 
accretion phase. 

The remainder of this paper is organized as follows. In 
Sect. |2 we explain the numerical method and physical assump- 
tions for constructing 2-D white dwarf models. Basic model 
properties are described in Sect. [3] and, in the next section 
(Sect. 0J, we present the quantitative criteria for thermonu- 
clear explosion and electron-capture (EC) induced collapse of 
white dwarfs in terms of the white dwarf mass and total an- 
gular momentum. In Sect. [5] the 1-D models of Paper I are 
compared to the 2-D models of the present study. We discuss 
implications of our results for the evolution of accreting white 
dwarfs, and suggest new evolutionary scenarios for SNe la 
progenitors in Sect [6] Possible consequences for the diversity 
of SN la are elaborated in Sect. [7] while implications of the 
EC induced collapse of differentially rotating white dwarfs are 
briefly discussed in Sect. [8] We summarize our main conclu- 
sions in Sect. [9] 

2. Constructing 2-D models 

2.1. Basic assumptions 

Rotating stars are called barotropic if the following two condi- 
tions are satisfied (Tassoul 1978): 

P = P(p) and ^=0, (1) 
oz 

where P is the pressure, p the density £2 the angular veloc- 
ity and z the distance to the equatorial plane, respectively. 
Degenerate cores of accreting white dwarfs are thought to 
be barotropic on two accounts. First, the equation of state 
of highly electron-degenerate matter is given by P — P(p). 
Second, perturbations which create + will decay on 
a short time scale in degenerate matter, as discussed in 
Kippenhahn & Mollenhof \ 19741 1 and Durisen ( I1977I I. Since 
the mass of the non-degenerate outer envelope in an accret- 
ing white dwarf is usually small (Paper I; Yoon et al. 2004 1, 
the barotropic condition may be a good approximation for de- 
scribing the whole white dwarf. For constructing axisymmetric 
2-D white dwarf models, we adopt cylindrical rotation as de- 
scribed below, and the equation of state of a completely degen- 
erate electron gas, as in Ostriker & Bodenheimer I1968> : 

P = A[x(2x 2 - 3)(x 2 + 1) 1/2 + 3 sinlT 1 x] ; p = Bx 3 , (2) 

where A = 6.01 xlO 22 dynes cirT 2 andB = 9.82xl0 5 // e g cirT 3 . 
Here we adopt fi e — 2 for all calculations. 

2.2. The rotation law 

According to Paper I, the rotation law in a white dwarf which 
accretes matter and angular momentum with accretion rates 
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higher than 10~ 7 M /yr is critically affected by the angu- 
lar momentum transport due to the dynamical shear instabil- 
ity (DSI) in the inner core, and due to the secular shear in- 
stability, Eddington-Sweet circulations and the GSF instabil- 
ity in the outer layers. In particular, the shear rate in the in- 
ner core remains near the threshold value for the onset of the 
DSI (:= cr D si,crit) during the accretion phase (see Paper I for 
more details). Another important feature of the rotation law is 
the presence of an absolute maximum in the angular velocity 
profile just above the shear-unstable core, which occurs due 
to the compression of the accreted layers with time (cf. also 
Durisen ( 1977 1 and see Fig. 9 in Paper I). From these con- 
straints, we can construct a realistic rotation law as follows. 

We define s p as the position of the maximum angular veloc- 
ity. As discussed in Paper I, above this position the buoyancy 
force is strong enough to assure stability against the DSI. In the 
white dwarf models in Paper I, this position is linked to layers 
with a density as low as several percents of the central density 
(p c ). Therefore, we determine s p such that 



p(s ■ 



^ = 0)=/„Pc 



(3) 



In our models, we use f p = 0.1 in most cases. Then, O(s) is 
determined in the range of < s < s p from the fact that, in the 
inner core, the shear rate is given by 



d£l(s) 
dins 



- /shC DSI.crit 



(4) 



Here, s is the distance from the rotation axis in the cylindri- 
cal coordinate and the dimensionless parameter / s h (< 1.0) de- 
scribes the deviation of the shear rate from cr D si jCr i t . We use 
/ s h = 1.0 except for some test cases. The threshold value of 
the shear rate for the onset of the DSI is given as cr D si jCr it = 
tJN 2 /Ri c , where N 2 is the Brunt-Vaisala frequency and R[ c the 
critical Richardson number (see Paper I for more details). From 
this we get 



Q(s) = £2 C + 



C S /shC DSI ,crit 

Jo 



ds' , 



(5) 



where Q c is the angular velocity at the center. In order to ensure 
that d£l/ds = at s = s p , which is the case in 1-D models of 
Paper I, we change the values of cr D si, c rit around s p such that 



C" DSI ; crit(s) - 



^/N T /R~ if < S < Sb 

y/N 2 /Ri c (s - s p )/(ib - s p ) if s b < s < s p 



,(6) 



where s b is a certain location near to s p . In this study, we define 
Sb as the fifth neighboring grid point from s p (while 149 grid 
points in ^-direction are used). The calculation of N 2 is car- 
ried out by the use of the equation of state given by Blinnikov 



et al. (1996 1, using the equatorial values of p and assuming 
T = 10 8 K which is typical in the inner core of the accret- 
ing white dwarf models in Paper I. For the critical Richardson 
number, R\ c = 0.25 is used. The integration in Eq. (|5} is car- 
ried out numerically by fitting / s hO"Dsi,crit/s with the Chebyshev 



polynomial (cf. Press et al. 1986). 

At the surface the white dwarf equator is assumed to rotate 
at a certain fraction (/k) of the Keplerian value, such that 



where R e is the equatorial radius of the white dwarf, and 
An accreting white dwarf may be ex- 
pected to rotate near the critical velocity at its surface, since ac- 
cretion through a Kelperian disk will carry angular momentum 
with the local Keplerian value. Therefore, we use /k = 0.95 in 
most cases. 

It remains to determine Q(s) between s p and R e . We impose 
an analytical form to Q(s) in this region as follows: 



£l(s)/Q K (s) = Q(s p )/Q K (sp) + C(s - s p ) a , 



(8) 



where Q.k(s) is the local Keplerian value. The constant C can 
be determined with a given value of a as 



C 



/ K - n(s p )/n K ( Sp ) 
(R e - s p y 



(9) 



The white dwarf structure is not significantly affected by the 
choice of the index a, because of the constraints imposed by 
£l(s p ) and Q(/? e ) at each edge. Furthermore, the density in the 
outer region (s > s p ) is lower than the inner core, and the role 
of the rotation law of this region for the white dwarf structure is 
less important than that of the inner core. In the present study, 
a - 1 .2 is usually employed. 

In summary, we determine the rotation law as follows. In 
the inner region (0 < s < s p ), Q.(s) is determined by the dy- 
namical shear instability ( Eqs. (0 and (|6}) for a given Q c . In 
the outer region (s p < s < R e ), Q.(s) is determined by Eq. 
where C is given by Eq. l|9}. 

2.3. SCF method 

For constructing 2-D white dwarf models with the prescribed 
rotation law, we use the so-called self-consistent field (SCF) 
method which has been widely employed for modeling rotating 
stars (Ostriker & Mark 1968 Hachisu 1986). This method can 
be most easily employed with barotropic stars, which satisfy 
the relation p = p( v F), where the total potential T* is the sum of 
the gravitational potential <t> and the centrifugal potential 0: 



r) = <D(p, r) + ®(s) , 



where 
<50u, r) - G 



J Ir-r'l 



and 



-f 

Jo 



ci(s'y s 'ds' 



(10) 



(11) 



(12) 



Q(fl e ) = / K Q K CR e ) , 



(7) 



Here, //(:= cos 9) and r represent the spherical coordinates. The 
distance from the rotation axis s is given by s = r(l - fi 2 ) 1 ^ 2 . 
In our study, the gravitational potential <£ is calculated using 
Simpson's formula following Hachisu ( 1986 1. For the centrifu- 
gal potential 0, the integration in Eq. (I12> is carried out by 
approximating the integrand by Chebyshev polynomials. 

The SCF method then proceeds as follows. Once we have 
T* from a given trial density distribution p(/i, r) and a given 
angular velocity profile, a better value for the density p can 
be given analytically as a function of *F (see Ostriker & 



4 



Yoon & Langer: Rapidly rotating massive white dwarfs 



Bodenheimer 1968 1. This new density distribution in turn gives 
a new value of x ¥. This procedure is iterated until a good con- 
vergence level is achieved. In the present study, the iteration 
continues until the following conditions are satisfied: 



vjjn 



■4" 



old 



< 5-10 , and 



C- < C > 



<C> 



< 5-l(T 5 ,(13) 



where C, which should be constant in ideal equilibrium struc- 
tures, is defined as 



C 



-I 



p" l dP - O - 



(14) 



and < C > denotes its mean value. One additional condition is 
imposed by the virial relation such that 



\2T- w + 3n\ 

w 



< 10" 



(15) 



where T and W denote rotational- and gravitational energy, re- 
spectively, and n the volume integral of the pressure. Note that 
here we define gravitational energy as a positive value for con- 
venience. Hereafter, we will refer to the above explained pro- 
cedure as SCF iteration. 

2.4. Further details 

We use the maximum density (p max ), the gravitational constant 
(G) and the equatorial radius at the surface of the white dwarf 
(R e ) for the basis of dimensional variables for numerical cal- 
culations, as in Hachisu ( 1986 1. With this choice, the physical 
variables can be expressed in a dimensionless form such that 

P = P/Pmax, S = S/R e , j = jl VOTw^e, Q = Q/ V^Pmax, 

*F = G/?eP m ax, and so on. As pointed out by Hachisu ( 1 986 1, 
using a fixed value of p max instead of the total mass M makes 
the numerical convergence faster, especially because the equa- 
torial radius R e can be given analytically as 



Re 



1/2 



(16) 



where indicates the boundary value. 

Hachisu also chose the axis ratio as an initially fixed param- 
eter for determining the equilibrium structure. This approach is 
possible only when the rotation law is given analytically, and 
can not be applied for the rotation law defined in Sect. 12.21 
Instead, we choose the central angular velocity Q c as an initial 
parameter. Therefore, in our models, initially fixed maximum 
density p max and central angular velocity £2 C determine the 
equilibrium structure of the white dwarf, together with given 
values of / p (Eq.0, / sh (Eq.0, / K (Eq.Q and a (Eq.|9J. More 
specifically, at a given p max , the total angular momentum J and 
the white dwarf mass M increase with increasing Q c . 

It is found that, when the centrifugal potential © is calcu- 
lated directly from Q(s) given by Eq. l|5} and (|SJl, models do not 
converge well if the ratio of the rotational energy to the gravita- 
tional energy (T/W) becomes rather high (> 0.07). This is due 
to the complexity of the prescribed rotation law. Usually, our 
numerical code works very efficiently even when T/W > 0.2 if 
the rotation law is a well defined function of position as in the 



case of v-constant rotation or j-constant rotation (see Eqs. El 
and^Jbelow). We thus pre-determine the specific angular mo- 
mentum j(s) from a trial density distribution before starting 
the SCF iteration such that j(s) = Cl{s)s 2 , where Q. is obtained 
from Eq. © and (|8j. We then introduce a new rotation law 
u>(s) such that a>(s) = j(s)/s 2 . Although Q. and a> are identical 
in the beginning of the SCF iteration, the two deviate each other 
during the SCF iteration, since Q. depends on both density and 
position (see Eq.[5]and|6jl while u>, given by the pre-determined 
j(s), depends only on position. Once the SCF integration is 
completed, we calculate again Q. by the use of Eq. and © 
from the converged model, and determine a new j(s). Then, the 
SCF iteration is carried out again with the new j{s). This pro- 
cedure is repeated until Q. and u converge to each other, within 
an accuracy of 3%. 

2.5. Model sequences 

Hereafter, we will refer to the rotation law given by Eq. $5^ and 
Eq (|8j to accreting-white-dwarf rotation (AWD rotation). For 
comparison, we employ rigid, v-constant and y'-constant rota- 
tion laws, which were also considered by Hachisu ( 1986 1. 

Model sequences which are constructed in the present study 
with different parameters and different rotation laws are sum- 
marized as follows (see also Tables 1 and 2, below): 

1 . AWD rotation: Q. given by Eq. l|5} and © 

- Sequence AA: / sh = 1.0, / K = 0.95, f p = 0.10, a = 1.2 

- Sequence AB: / sh = 1.0, f K = 0.95, f p = 0.05, a = 1.2 

- Sequence AC: / sh = 1.0, f K = 0.95, / p = 0.10, a = 1.0 

- Sequence AD: / sh = 1.0, / K = 0.95, / p = 0.20, a = 1.2 

- Sequence AE: / sh = 0.5, f K = 0.50, f p = 0.10, a = 1.2 

- Sequence AF: / sh = 0.1, f K = 0.95, f p = 0.10, a = 1.2 

- Sequence AG: / sh = -1.0, f K = 0.6, / p = 0.01, a = 0.9 

2. rigid rotation (Sequence R): Q = constant 

3. y'-constant rotation (Sequence J): 



Q = ;o/(0.01fl e 2 + s 2 ) . 
4. v-constant rotation (Sequence V): 
Q = v /(0.01/? 2 + i 2 ) 1/2 . 



(17) 



(18) 



Here, sequences AA, AB & AC are intended for describing 
the properties of 1-D models in Paper I. Sequences AD, AE & 
AF are to consider uncertainties due to changes of the given pa- 
rameters (see below). Sequence AG is to investigate the white 
dwarf structure when the inner core rotates at the shear rate of 

-0"DSI ; crit- 

We consider model sequences at 7 different values of p max : 
1.0, 2.0, 5.0, 8.0, 20, 40.0, and 100 xlO 8 gem- 3 . Every se- 
quence is given a subindex (a, b, c, d, e, f and g), in order 
to indicate the corresponding p max . For instance, AAa denotes 
the sequence with / S h = 1.0, /k = 0.95, f p = 0.10, a = 1.2, 
and with p max = 1.0 x 10 8 g cirT 3 , etc. Note that in all mod- 
els with the AWD- or rigid rotation law the maximum den- 
sity p max corresponds to the central density p c . However, they 
may differ from each other for v-constant and /'-constant rota- 
tion laws when the ratio of the rotational energy to the gravi- 
tational energy (T/W) becomes sufficiently high. Within each 
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Fig. 1. Upper panels: Angular velocity as a function of the mass coordinate in white dwarf models (AAal6, AAa24 & AAa31 
in the left panel, and AAelO, AAe24 & AAe37 in the right panel). Lower panels: Angular velocity normalized to the local 
Keplerian value as a function of the mass coordinate, in the corresponding models. 



sequence, white dwarf models with different values of Q. c are 
constructed. Since the mass budget in binary systems which 
are considered for the Single Degenerate scenario for SNe la 
is limited, the maximum possible mass that a white dwarf 
can achieve by mass accretion could not significantly exceed 
2.0 M (Langer et al. 2000l. Therefore, we do not construct 
models with M > 2.1 M Q . 

3. Basic model properties 

Properties of selected models with the AWD rotation law are 
presented in Tables ^ ar, d ID As mentioned earlier, Q pro- 
files in the models with / s h = 1.0, / p = 0.05 or 0.1, and 
/ K = 0.95 (Seqs. AA, AB, & AC) resemble closest the 1-D 
models in Paper I. The models with p c = 2 x 10 9 g cm" 3 repre- 
sent the pre-explosion stages, since carbon burning is expected 
to occur at roughly this density in accreting white dwarfs with 
M > 10" 7 M yr 1 (e.g. Nomoto fT982l 

Fig-fflshcws angular velocity profiles of six different mod- 
els, as indicated. Here, the mass coordinate m(s) is defined such 
that m(s) - 2n f. f ^pis' ,z)dzds' . These profiles have simi- 
lar shapes to those in 1-D white models in Paper I. A detailed 
comparison to 1-D models is performed in Sect. [5] Note that 
models with p c = 2 x 10 9 g cirT 3 have higher values of Q. than 
those with p c = 1 x 10 8 g cirT 3 at a given mass, due to their 



stronger compactness. This tendency is reflected in the values 
of the moment-of-inertia weighted mean of Q in Tabled 

Iso-density contours of the corresponding models with p c = 
2 x 10 9 g cm" 3 are given in Fig. [2 As the value of T/W in- 
creases, the white dwarf structure deviates from the spheroidal 
shape more and more, and the ratio of the polar radius to 
the equatorial radius decreases. However, the inner core with 
p > 0.05p c , which contains more than 75% of the total mass, 
remains fairly spheroidal even when T/W = 0.1 1. 

Even though our 2-D models are designed to resemble 
the results from evolutionary calculations, the parameter range 
which describes them is finite due to two effects. Firstly, vari- 
ations in the evolution itself, i.e. of the accretion rate in bi- 
nary systems, or of the initial white dwarf mass, give rise to 
variations, say, in the pre-supernova modes (i.e., models with 
p c = 2 x 10 9 g, cm" 3 ); i.e., the value of / p has been shown to de- 
pend on both in Paper I. Secondly, some physical processes de- 
scribed in the evolutionary models are uncertain. For instance, 
the exact value of the critical Richardson number is still uncer- 
tain, and instead of R hC = 0.25 that is used in the present study, 
^i,c = 1-0 is also often adopted in the literature (e.g. Hirschi et 
al. 2004), with which the shear rate in the inner core of white 
dwarfs would become weaker than assumed in sequence AA. 
Therefore, we constructed 2-D models in which we varied all 
relevant parameters (Table 
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Fig. 2. Iso-density contours of white dwarf models AAel4, 
AAe24 and AAe37. These models have the same central den- 
sity of p c — 2 x 10 9 g cm 4 , but different masses and angu- 
lar momenta. The contour lines denote levels of 0.8p c , 0.6p c , 
0.4p c , 0.2p c , O.lpc, 0.05p c , l(T 2 p c , 10 4 p c , 10" 4 Pc and l(T 5 p c , 
from the central region to the outer layers. 
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Fig. 3. Upper panel: Angular velocity as a function of the 
mass coordinate in white dwarf models AAe20 (solid line), 
ABe20 (dotted line), ACel9 (short-dashed), AFe8 (dashed- 
three-dotted line) and AGe3 (long-dashed line). Lower panel: 
Specific angular momentum in the corresponding models of the 
upper panel. 



Fig. [3] shows Q-profiles in 5 different models (AAe20, 
ABe23, ACel9, AFe8 & AG3), which have the same central 
density (p c = 2 x 10 9 g, cm 4 ) and similar masses (~ 1.58 M Q ), 
but different sets of f p , / s h, /k and a (see Tables [2 and A 
comparison of model AAe20 with model ABe23 shows the 
influence of / p . For smaller / p , the Q-peak is located further 
outward. As a result, the inner core rotates more slowly and 
the outer region more rapidly, for fixed mass and central den- 
sity. The influence of a can be understood by comparing model 
AAe20 with model ACel9: the surface region rotates more 
rapidly and the inner core more slowly in model ACel9. This 
is because, with a smaller a, the role of the centrifugal force 
becomes stronger in the region of s > s p (see Eq. [8jl. Model 
AFe8 has a weaker shear rate in the inner core than model 
AAe20. As a consequence, this model needs a higher value of 
Q c than model AAe20, to retain the same central density at a 
given mass. In model AGe3, where the Q-gradient is negative 
throughout the white dwarf, the inner core rotates much faster 



than the other models. In Fig.0] models with /'-constant and 
v-constant rotation laws are compared with models AAe20 and 
AGe3. Note that differential rotation in those models is much 
stronger than in model AGe30, which means that the shear rates 
for such exotic rotation laws are well above the critical value 
for the onset of the dynamical shear instability, cr D si jCr i t . 

The specific angular momentum profiles in Fig.[3]show an 
interesting feature. If a white dwarf model has lower values 
of j in the inner core than the other models do, it has higher 
values of j in the outer layers. As remarkable consequence, 
all the white dwarf models displayed in Fig. [3] have a similar 
amount of total angular momentum (see Table^andlJ}, despite 
their different angular momentum distributions. Fig [5] where 
the total angular momentum J of models in different model se- 
quences (AAe, ABe, ACe, AFe, AGe & Re) is plotted in the 
M — J plane, indeed shows that the total angular momentum, 
for a given central density, is a function of the white dwarf mass 
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Fig. 5. Total angular momentum of white dwarf models as a 
function of the total mass of the white dwarf, in sequences 
AAe, ABe, ACe, AFe and AGe. 



only, regardless of various possible rotation profiles. Although 
the scatter in the value of J becomes larger as the white dwarf 
mass increases, it remains within 10 % from the mean value 
even when M > 1.9 M . Even the rigidly rotating white dwarf 
models (M < 1 .47) follow the general trend very well. We con- 
clude that the 7-M-law predicted by our models is not sensitive 
to evolutionary details or uncertain physical processes. 

For comparison, we also give results with j'-constant and 
v-constant rotation laws in the same M - J plane. Note that 
models with these exotic rotation laws follow very different re- 
lations. However, as mentioned above, these models are shear 
unstable, and their steep angular velocity gradients should de- 
cay on a short time scale (see Paper I). Therefore, we can con- 
clude that M— J relation is not significantly affected by different 
O-profiles in general. We find that this is also the case for the 
dependence of other physical quantities on the total mass, such 



as for T/W and and for the binding energy, as discussed in the 
following section. 

4. Functional relations of physical quantities 

Here we derive analytical expressions for various physical 
quantities as functions of the white dwarf mass (M) and the 
central density (p c ), from models with the AWD-rotation and 
rigid rotation laws. 

4.1. M - J relation 

Fig-Elshows the total angular momentum of white dwarf mod- 
els with the AWD or rigid rotation law for 7 different cen- 
tral densities, in the M - J plane. White dwarfs with lower 
central densities have systematically larger angular momenta 
at a given mass, as expected from the role of the centrifugal 
force. We find that the relation between the total angular mo- 
mentum and the total mass of the white dwarf can be given 
in an analytical form, in the considered range of the central 
density (10 8 < p c (g cm 4 ) < 10 10 ) and the white dwarf mass 
(M < 2.0 Mq), as 

J(M;p c ) = d(p c ){l -exp (-0.2 [M-M NR (p c )f 48 )} 

+C 2 (Pc)[M-M NR (p c )] L3 [l0 50 ergs], (19) 

where 

Ci(p c ) = 20.800370 - 1 .5856256 log I0 p c , (20) 
and 

C 2 (p c ) = 11.397626- 0.97306637 log 10 p c . (21) 

Here, Mnr(p c ) denotes the mass of the non-rotating white 
dwarf at a given central density, and both M and Mnr are given 
in the units of solar masses. From non-rotating white dwarf 
models, we derive Mnr as a function of p c as 

M NR = 1.436 x 

[l - exp {-0.01316(log loPc ) 2 - 706 +0.2493 log I0 p c }] 
(forp c > 10 7 gcnT 3 ) . (22) 

We are particularly interested in the cases of p c = 2 x 
10 9 g cm 4 , p c = 4 x 10 9 g cm 4 and 10 10 g cm 4 , as they 
represent critical densities for the thermonuclear explosion of 
CO white dwarfs, and the electron-capture induced collapse 
of ONeMg and CO white dwarfs, respectively (see Nomoto & 
Kondo 199TJ. Hereafter, we refer to the total angular momen- 
tum at these densities as /sNia, •/ec.onm, and /ec,co, respec- 
tively: 



^SNia = J(M;p c = 2 x 10 9 g cm 3 ) , 
/ec,onm = J(M;p c = 4 x 10 9 g cm 4 ) , 
and 

/ec,co = J(M; Pc = 10 10 gcm 4 ). 



(23) 
(24) 

(25) 
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Fig. 6. Upper Panel: Total angular momentum of white dwarf 
models with the AWD rotation and rigid rotation laws, as a 
function of the white dwarf mass. Models at different cen- 
tral densities are plotted with different symbols as indicated in 
the figure. The solid lines denote the fitting function given by 
Eq. ( I19> . The dashed line gives the M - J relation of rigidly 
rotating white dwarfs at critical rotation. Lower Panel: Specific 
angular momentum of white dwarf models with the AWD rota- 
tion and rigid rotation laws, as a function of the white dwarf 
mass. Models with p c = 2 x 10 9 g cirT 3 , 4 x 10 9 g cm" 3 
and 10 10 g cirT 3 are plotted with different symbols as indi- 
cated in the figure. The solid lines represent J(M; p c )/M, where 
J(M; p c ) is given by Eq. dl9l . The dashed line gives the M - j 
relation of rigidly rotating white dwarfs at critical rotation. 



If a CO white dwarf has a larger amount of angular momentum 
than /sNia, a SN la explosion can not be induced. Similarly, 
an ONeMg white dwarf with / > /ec.onm will not experi- 
ence electron-capture induced collapse. On the other hand, a 
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Fig. 7. The critical angular momentum for thermonuclear ex- 
plosion of CO white dwarfs (/sNia; solid line), and electron- 
capture induced collapse of CO white dwarfs (7ec,co; dashed- 
dotted line) and ONeMg white dwarfs (7ec,onm; dashed line), 
as a function of the white dwarf mass. SNIa explosion is ex- 
pected when 7ec,co ~ J ~ Anis (hatched region). Electron- 
capture induced collapse is supposed to occur when J < 7ec,co 
for CO white dwarfs, and J < /EC.ONeMg for ONeMg white 
dwarfs, respectively. 

SN la explosion is expected when 7ec,co ~ J ~ -/sNia- If car- 
bon ignites at the center when J < /ec.co, the outcome will 
be electron capture induced collapse. These criteria for SN la 
explosion and EC induced collapse are summarized in Fig. [71 

Alternatively, the critical CO white dwarf mass for a SN la 
explosion (MsNia) an d ONeMg or CO white dwarf mass for 
EC induced collapse (Mec.onm or Mec.coX can be given as a 
function of /, for M < 2.0 M , as 



MsNia = 1-376 - 0.024/ - 83 + 1 - exp(-0330J i : ') , 



,1.5% 



(26) 



Meconm = 1-395 ■ 



0.040/ - 83 + 1 



• exp(-0.4167 L5 ) , (27) 



and 

-Meqco 



1.412 



0.0427 83 + 1 - exp(-0.5197 15 ) 



(28) 



I.e., exploding/collapsing white dwarfs are expected to have 
roughly these critical masses for a given total amount of an- 
gular momentum. 

4.2. M-T/W relation 

The ratio of the rotational energy to the gravitational energy 
(T/W) is plotted as a function of the white dwarf mass in Fig. [8] 
for models with the AWD rotation and rigid rotation laws. For 
a fixed central density p c , the values of T/W follow very well 
a one-dimensional relation; the scatter is much smaller than in 
the corresponding diagram showing the total angular momen- 
tum (Fig. [6}. This relation can be given as a function of M and 
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Fig. 8. Ratio of the rotational energy to the gravitational energy 
(T/W) in white dwarf models with the AWD rotation and rigid 
rotation laws, as a function of the white dwarf mass. Models at 
different central densities are plotted with different symbols as 
indicated in the figure. The solid lines denote the fitting func- 
tion given by Eq. i29\ . 

Pc, as 

T/W(M;p c ) = 0.1{l-exp(-C 3 (p c )[M-M NR (p c )])} 

+C 4 (p c ) [M - M NR (p c )] , (29) 

where 

C 3 (p c ) = 407.71741 - 168.99280 log 10 p c 

+ 26.449872(log 10 p c ) 2 - 1.8421546(log 10 p c ) 3 (30) 



+ 0.048 12 1239(log 10 p c ) 4 



and 

C 4 (p c ) 



-3.2305060+ 1.5464029 log 10 p c 
0.26049998(log 10 p c ) 2 (31) 
0.01881 1083(log 10 p c ) 3 - 0.00049405568(log I0 p c ) 4 



Note that, in the considered mass range (M < 2.0 M ), only the 
models with p c = 1 x 10 8 g cirT 3 and 2 x 10 8 g cm 4 reach the 
canonical value of 0. 14 for the onset of the bar-mode instability 
(e.g. Durisen & Imamura ll977> . 

4.3. Binding energy 

The binding energy (BE :- W - U - T) is essential for the 
energy budget of exploding or collapsing white dwarfs. Fig.|5] 
shows binding energy and specific binding energy of our white 
dwarf models with p c = 2 x 10 9 g cm 4 and p c = 10 10 g cm 4 , 
which may represent the pre-explosion and pre-collapse stages 
of CO white dwarfs, respectively, and with 4 x 10 9 g crrT 3 , 
which refers to the pre-collapse stage of ONeMg white dwarfs. 
Both, binding energy and specific binding energy are monoton- 
ically increasing functions of the white dwarf mass. 
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Fig. 9. Upper panel: Binding energy of white dwarf mod- 
els with the AWD rotation and rigid rotation laws for p c = 
2xl0 9 gcm 3 4x 10 9 g cm 3 , and 10 10 g cm 3 , as function of 
the white dwarf mass. The solid lines denote the fitting func- 
tion given by Eq. i32\ . Lower panel: Specific binding energy of 
white dwarf models with the AWD rotation and rigid rotation 
laws for p c = 2 x 10 9 g cm 4 , 4 x 10 9 g cm 4 , and 10 10 g cm 4 . 
The solid lines are from the fitting function given by Eq. d32i . 
divided by M. 

We constructed a fitting formula for the binding energy as 
a function of the white dwarf mass, which is applicable in the 
considered range of central densities (1 x 10 8 < p c [g cm 4 ] < 
10 10 ): 

BE(M;p c ) = BE NR (p c ) 

+C 5 (p c )[M-M NR (p c )] 103 [l0 50 erg], (32) 

where 

C 5 (p c ) = -370.73052+ 132.97204 log 10 p c 

- 16.11 703 l(log 10 p c ) 2 + 0.66986678(log 1() p c ) 3 . (33) 

The binding energy of non-rotating white dwarfs BE^r is given 
by 

BE NR (p c ) = -32.759747 + 6.7179802 log 10 p c 
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Fig. 10. Upper panel: Ratio of rotational to gravitational en- 
ergy (T/W) as a function of the white dwarf mass. The thin 
solid line corresponds to the M - T/W relation of rigidly ro- 
tating 2-D white dwarf models at critical rotation. Thick solid 
lines are evolutionary tracks of the 1-D sequences C2 and C10 
in Paper I. Labels at filled circles on those tracks denote cen- 
tral densities in the unit of 10 8 g cirT 3 . The M - T/W re- 
lation given by Eq. J29l > is also plotted for different densi- 
ties: p c /[10 8 g cm" 3 ] = 1.0 (dotted line), 2.0 (dashed line), 4.0 
(dashed-dotted line), 8.0 (dashed-three-dotted line), 20.0 (thick 
dashed line) and 40.0 (long-dashed line). Lower panel: Same 
as in the upper panel, but in the M - J plane. 



-0.287 17609(log 10 p c ) 2 



[l0 50 erg] 



(34) 



5. Comparison with 1-D models 



In Fig. we compare our 2-D results with with two time- 
dependent 1-D model sequences of Paper I in the M—T/W and 
M— J planes. These sequences simulate white dwarfs which ac- 



crete with constant accretion rates of 5 x 10~ 7 M yr _1 , for an 
initial white dwarf mass of 0.8 M (sequence C2) and 1.0 M 
(sequence C10). Fi sure ITOl shows that the models of sequence 
C10 have similar T/W values at a given mass and a central den- 
sity as the 2-D models as long as M < 1.6. However, they devi- 
ate significantly for higher masses. The models of sequence C2, 
which have a higher angular momentum at a given mass com- 
pared to those of sequence C10, deviate from the 2-D models 
already when M > 1.4. We conclude that the 1-D models can- 
not be regarded as accurate when T/W > 0.06, which occurs 
after accreting about 0.4 M in both cases. 

The inconsistency between 1-D and 2-D models at T/W > 
0.06 is not surprising, because in 1-D models, the effect of the 
centrifugal force on the white dwarf structure cannot be accu- 
rately described when a significant fraction of the white dwarf 
rotates at above ~60% of the Keplerian value (cf. Paper I). In 
particular, the white dwarf radius R is underestimated in this 
situation. In the 1-D models, a certain amount of angular mo- 
mentum is assumed to be accreted onto the white dwarf, and the 
consequent angular velocity is determined as Q ~ J /(MR 2 ). 
From T ~ MQrR 2 ~ J 2 /(MR 2 ) and W ~ M 2 /R, we have 
T/W ~ J 2 /(M 3 R). Therefore, at given J and M, an underesti- 
mate of R leads to an overestimate of T/W. 

It is remarkable that 1-D and 2-D models show an extraor- 
dinary agreement in the M - J plane throughout the consid- 
ered mass range (Fig. [K)] lower panel). However, the impli- 
cation is not that 1-D models can describe the white dwarf 
structure accurately. It should be rather considered as coin- 
cidence due to the following reason. The underestimate of 
R in 1-D models leads to an overestimate of p c , since it is 
p c ~ M/R 3 . On the other hand, the central density should be 
reduced due to the overestimate of the centrifugal force, as the 
ratio of the centrifugal force to the gravitational force is given 
by FJF g ~ D?R/(GM/R 2 ) ~ J 2 /(M 3 R). These two opposite 
effects might just compensate each other. 

In any case, the consistency between 1-D and 2-D models 
in the M - J plane leads us to retain one of the most impor- 
tant conclusions of Paper I: If white dwarfs accrete mass and 
angular momentum efficiently, they might hardly reach central 
carbon ignition during the accretion phase. 

6. Implications for the SNe la progenitor evolution 

Our 2-D models of rotating white dwarfs show that various 
physical quantities can be given as a function of the white 
dwarf mass, as long as the shear rate inside the white dwarf 
is not stronger than the critical value for the onset of the 
DSL In particular, the M - J relations of white dwarfs at 
p c = 2 x 10 9 g cm" 3 (/sNia) and p c = 10 1() g cm" 3 (7 E c,co) 
provides a quantitative criterion for the supernova explosion 
of rotating white dwarfs, as summarized in Fig. [7] Here, we 
discuss possible evolutionary paths of accreting white dwarfs 
towards a supernova explosion. 



6.1. The angular momentum problem 

If a non-rotating white dwarf with the initial mass of Mi n ; t ac- 
cretes mass with a specific angular angular momentum of the 
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Keplerian value (i.e., = y/GMR), the total amount of ac- 
creted angular momentum when the white dwarf reaches mass 
M f is 

J = j jdM ^ VGA [m] 12 - M^] (35) 

,4xl0-(^) 1/2 f(^f-(^n (ergs, 
\0mR Q J [\M Q J \ Ma ) J S 

In binary systems which are considered for the so-called Single 
Degenerate scenario for SNe la, the mass budget is limited, 
and the maximum possible mass which a CO white dwarf 
can achieve by mass accretion is about 2.0 M (Langer et 
al. |3T)00l. Therefore, white dwarfs which are potentially able 
to perform a SN la explosion will have masses in the range of 
1.4 < M/M Q < 2.0. Given that initial masses of 0.8... 1.0 M are 
usually considered in this context (e.g. Langer et al 2000 1, the 
value of J according to the above equation will reach 2.6...8.1 x 
10 50 ergs when the white dwarf mass reaches 1.4. ..2.0 M . 
This is certainly much larger than /sNia (see Fig.0, and cen- 
tral carbon ignition can not occur in white dwarfs with such 
large angular momenta. 

For accreting white dwarfs to end in a SN la, therefore, 
either smaller amounts of angular momentum need to be ac- 
creted, and/or a large fraction of the accreted angular momen- 
tum should be lost again. This means that the detailed evolu- 
tionary paths of SNe la progenitors will depend on the effi- 
ciency of the angular momentum gain and on the time scale of 
the angular momentum loss. 

6.2. Efficiency of the angular momentum gain 

In Paper I and also in Uenishi et al. ( 2003 ) and Saio & Nomoto 
( 2004 1, the outward angular momentum transport from the crit- 
ically rotating star into the accretion disk is considered to al- 
leviate the angular momentum problem, based on the work by 
Paczyhski ( 1991 1 and Popham & Narayan ( 1991 1. For instance, 
in Paper I, if the equatorial surface of a white dwarf model ro- 
tated at the Keplerian value, angular momentum was assumed 
to be transported outward from the white dwarf into the ac- 
cretion disk, while otherwise angular momentum of a certain 
fraction of the Keplerian value was allowed to be accreted. 
However, the results of Paper I showed that even with this limi- 
tation, the efficiency of the angular momentum gain (hereafter, 
7 acc -efficiency) is usually too high for producing SNe la dur- 
ing the accretion phase: only when the amount of the accreted 
angular momentum was limited to 30% of the Keplerian value, 
white dwarfs with a high initial mass (M; n j t = 1.0 M ) could 
reach carbon ignition at the center, which occurred at about 
M ^ 1.8 M Q (sequences C10, CI 1 & C12 in Paper I; see also 
Fig. [I0l. 

It should be noted that the angular momentum accretion 
efficiency might be below 100% on two accounts. Firstly, 
the models of Paper I neglected thermonuclear shell burning. 
Eddington-Sweet circulations will be induced in the presence 
of hydrogen and helium shell burning, which will accelerate 
the outward transport of angular momentum in the white dwarf 
envelope where the Q gradient is negative. Accordingly, the 



surface of the white dwarf will rarely deviate from rotating crit- 
ically, and the angular momentum gain from the accretion disk 
will be thus more severely restricted. 

Secondly, as pointed out by Saio & Nomoto (20041, the 
models of Paper I applied outward angular momentum trans- 
port into the accretion disk only for limiting the angular mo- 
mentum gain when the white dwarf surface reaches critical ro- 
tation, but did not allow for the possibility of removing angu- 
lar momentum from the white dwarf into the accretion disk. 
Simulations in Saio & Nomoto (2004i who did consider the 
latter possibility show that the total angular momentum in the 
accreting white dwarf can not only increase, but also decrease, 
depending on the change of the surface conditions. Systematic 
studies of the angular momentum accretion efficiency are there- 
fore highly desirable to improve our understanding of the evo- 
lution of SN la progenitors. 

6.3. Gravitational wave radiation 

As discussed in Paper I, rotating white dwarfs can be suscepti- 
ble to the bar-mode or to the r-mode instability, which will lead 
to loss of angular momentum via gravitational wave radiation. 
Comparison of Figs.|6]and [^indicates that J > 4 x 10 5() erg s 
is necessary for the white dwarf to reach T/W = 0.14, which 
is the canonical value for the onset of the bar-mode instability. 
Although Paper I described 1-D models with higher values of 
T/W, this is due to the overestimate of T/W in the 1-D approx- 
imation of the stellar structure, as discussed in Sect. [5] In terms 
of angular momentum, the 1-D models in Paper I never accrete 
more than 2.5 x 10 50 erg s, implying that T/W = 0.14 might 
hardly be achievable. We conclude that the bar-mode instabil- 
ity may be irrelevant on the SN la context, and only the r-mode 
instability may be relevant in accreting white dwarfs. 

Lindblom (1999 1 derived the time scale for the growth 
of r-mode instability in rigidly rotating stars as l/r r = 
K j p(r£l) 6 dr, where K = (2^/25)(4/3) 8 G/c 7 . Using the mean 
values for density, radius and angular velocity (p, R and £2), this 
time scale can be given as l/r r * KpilR 1 /7 = (3K/4tt)MD. 6 R 4 . 
Since we have Q. = J/(kMR 2 ), T = (k/2)MQ. 2 R 2 and W = 
yGM 2 /R, where k and y vary according to the stellar structure, 
the time scale r r can be expressed in terms of M, J and T/W as 
the following: 

log 10 Tr[yr] = -3.94+101og 1() (//10 50 )- 191og 10 (M/M G ) 

-8 log(77W) - 8 log 10 r -2 logic * ■ (36) 

In our white dwarf models, k and y are in the range of 0. 1 ~ 
0.14 and 1.0 ~ 1.6, respectively. 

In Fig. ^2 li nes of constant r r are plotted in the M— J plane, 
in the range of 1.4 < M/M Q < 1.9 and 0.5 < 7/[10 50 erg s] < 
2.0. Here, the M - J and M - T/W relations given by Eq. ( fT9l 
and J29i are used, with the fixed values of y — 1.0 and k — 0. 14. 
Values of r r are found to be about ~ 10 3 ... ~ 10 5 yr in the con- 
sidered range of M and J. This estimate indicates that the r- 
mode instability might have enough time to operate in white 
dwarfs during the accretion phase, if M =~ 10~ 7 M yr _1 , 
which is typical for the Single Degenerate scenario for SN la 
progenitors. 
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Fig. 11. Schematic supernova progenitor evolution scenarios in the M - J plane. Thick lines starting at a white dwarf mass of 
1 M Q with J = indicate various possible evolutionary paths of accreting white dwarfs: with angular momentum accretion 
and simultaneous angular momentum loss (A), without angular momentum loss during the accretion phase (B), with inefficient 
angular momentum accretion (C), and with maximum internal angular momentum transport efficiency (D); (see text for details). 
The region where 7ec,co ~ J ^ Jsnih is hatched, and the region where J < 7ec,co is colored. For M =s 10~ 7 . . . 10~ 6 M yr _1 , a 
supernova explosion is expected at the upper borderline of the hatched area, J = Jsma\ cf. Eq. t23i . The thin dashed line gives 
the M - J relation of rigidly rotating white dwarfs at critical rotation. Dotted lines denote lines of constant growth time of the 
r-mode instability (r r ; Eq.l36> and are labeled with the logarithm of the growth time in years. 



6.4. Evolutionary scenarios 

Let us define the accretion time scale as T acc , and the time scale 
for the loss of angular momentum due to r-mode instability as 
rj, respectively. 

6.4.1 . Case A: t t < r acc 

When a white dwarf accretes sufficient mass and angular mo- 
mentum, the r-mode instability begins to operate, and as a con- 
sequence the white dwarf loses angular momentum. Here we 
assume that, even though the angular momentum loss rate is 
small initially, it will increase during the white dwarf spin-up, 
to an extent that rj < r acc , i.e, that the loss of angular momen- 
tum occurs more rapidly than its accretion. Such a scenario is 
motivated by the r-mode angular momentum loss time scale as 
given by Eq.l36l 



In this case, the growth of the white dwarf angular momen- 
tum — and thereby the growth of its mass — is limited by the 
angular momentum loss; it ignites carbon at the center when J 
approaches /sNia, as illustrated by Path Al in Fig. \H\ The end- 
ing point of Path Al, which defines the explosion mass Mai, 
depends on both, the angular momentum loss time scale r r , 
and the angular momentum accretion efficiency. One dimen- 
sional evolutionary calculations with an angular momentum 
accretion efficiency as in Paper I show that rj r r results in 
M A i - 1-6 M and 1.5 M Q , for the initial masses of 0.8 M 
and 1.0 M (Knaap, Yoon & Langer, in preparation). 

In some binary systems, the mass budget may be limited 
and the mass transfer stops before the white dwarf reaches Mai , 
but at a mass higher than the classical Chandrasekhar mass. 
Then, the white dwarf evolves vertically downward in the M - 
J plane as illustrated by path A2. Since tj < 10 6 yr in this 
case, the central temperature at the moment when the white 
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dwarf reaches /sNia is high enough ( T 10 K) to ignite carbon 
at p c a 2 x 10 9 g cm -3 (see Paper I). Therefore, both Path 
Al and Path A2 lead to the production of SNe la. The masses 
of exploding white dwarfs through Path A are in the range of 
1.38M < M < Mai. Electron-capture induced collapse can 
not occur in Case A. 

6.4.2. Case B: tj > r acc 

If rj > r acc , the evolutionary path of the accreting white dwarf 
will critically depend on the angular momentum accretion ef- 
ficiency. If it is so high that the white dwarf angular momen- 
tum always exceeds /sNia during the accretion phase, the white 
dwarf mass will continue to grow until the mass transfer comes 
to an end. Once the mass transfer stops, carbon ignition will be 
delayed until the white dwarf loses enough angular momentum 
(Path B in Fig. II It . Its final fate will depend on the angular 
momentum loss time scale. If it is larger than the white dwarf 
cooling time (i.e., tj > 10 9 yr), carbon ignition will occur after 
the white dwarf core has crystallized — i.e., when J 7ec,co 
— , which leads to electron-capture induced collapse and re- 
sults in a neutron star (Nomoto & Kondo 1991 1. This case is 
illustrated by Path B2. 

If tj is short enough for carbon to ignite in the liquid state 
in the white dwarf core, a SN la explosion will occur (Path Bl). 
The mass of exploding white dwarfs can vary from 1.38 M Q to 
2.0 M , depending on the mass budget in the binary system. 
However, in this case, the large angular momentum loss time 
scale of more than 10 6 yr leads to an ignition density (pi gn ) sig- 
nificantly above ~ 2 x 10 9 g cm 4 . As SNe la at such high 
densities produce highly non-solar neutron-rich iron group el- 
ement ratios (Iwamoto et al. 1999), this scenario can not be 
considered as the major evolutionary path. 

6.4.3. Case C: tj > r acc and low 7 acc -efficiency 

If the angular momentum accretion efficiency is sufficiently 
low, white dwarfs can reach /sNia to produce a SN la dur- 
ing the accretion phase even if tj > r acc , as long as enough 
mass is supplied by the binary companion, as illustrated by 
Path CI. If the mass transfer stops when M > 1.38 M but 
before reaching /sNia, carbon ignition will be delayed, and its 
outcome will depend on the time scale for angular momentum 
loss, as in the Case B. If tj is shorter than the cooling time 
scale, a SN la explosion is expected soon after J reaches /sNia 
(Path C2). Otherwise, carbon burning will occur only when 
J =i 7ec,co and result in electron-capture induced collapse to 
form a neutron star (Path C3). The masses of exploding white 
dwarfs through Path C are between 1.38 M and the ending 
point of Path CI (M Cl ). 

6.4.4. Case D: Extreme J-transport efficiency 

Although Paper I showed that non-magnetic accreting white 
dwarfs are expected to rotate differentially, the existence of 
rigidly rotating white dwarfs, through extremely efficient trans- 
port of angular momentum by magnetic torques, can not be 



excluded. If this is the case, the white dwarf surface begins 
to rotate critically soon after the onset of mass accretion, as 
shown by the previous work (Piersanti et al. l2003al Uenishi et 
al. l2003l Saio & Nomoto 12004b 

The angular momentum of rigidly rotating white dwarf at 
critical rotation has an absolute maximum at about 1.13 M , 
and decreases as the mass increases further (see Fig. I10i . 
Therefore, once a rigidly rotating white dwarf reaches criti- 
cal rotation, it must lose angular momentum during the accre- 
tion when M > 1.13 M . The r-mode instability could remove 
angular momentum, if tj < r acc . However, in our rigidly ro- 
tating white dwarf models at critical rotation, the estimated 
growth time scale for the r-modes, r r , becomes smaller than 
10 6 yr only when M > 1.25 M . It is therefore question- 
able whether the r-mode instability can be an efficient mech- 
anism to remove angular momentum during the evolution of 
rigidly rotating white dwarfs. As mentioned earlier, also the 
bar-mode instability can not operate as the high required value 
of T/W ^ 0.14 can not be reached in rigidly rotating white 
dwarfs. Uenishi et al. {2003) and Saio & Nomoto (2004) con- 
sidered the viscous transport of angular momentum from the 
white dwarf into the accretion disk, so that a rigidly rotating 
white dwarf with M > 1.13 M may keep rotating critically 
with continuous accretion of mass, while angular momentum 
decreases progressively. 

Within this picture, the evolution of rigidly rotating white 
dwarfs will follow Path D in Fig. \H\ If the binary companion 
supplies sufficient mass to reach M ^ 1.46 M , a SN la ex- 
plosion occurs during the accretion phase (Path Dl). If mass 
transfer ceases before reaching M ^ 1.46 M but after M > 
1.38 M , carbon ignition at the center will be delayed until 
enough angular momentum is lost (Paths D2 and D3), with 
the fate of the white dwarf depending on Tj, as in the case 
B. Exploding white dwarfs through Path D have masses in the 
range 1.38 < M/M < 1.46 (Uenishi et al.. l2003l . 

6.4.5. Diversity in M S Nia 

The scenarios outlined above show that the masses of explod- 
ing white dwarfs (MsNia) may not be uniform, but are ex- 
pected to be diverse, with an explosion mass range depend- 
ing on the evolutionary path realized by nature. The minimum 
mass is the canonical value of 1.38 M , while the maximum 
possible mass (M max ) remains uncertain; Paths B and D re- 
sult in the largest (1.38 < M S Nia $ 2.0 M ) and smallest 
(1.38 < MsNia ~ 1-46 M ) diversity in mass, respectively. 
Paths A and C may lead to M max < 2.0 M , with the exact up- 
per mass limit depending on the angular momentum loss time 
scale and the angular momentum accretion efficiency. 

7. Observable consequences of 

super-Chandrasekhar mass explosions 

Recent results of 3-D numerical simulations of SN la explo- 
sion indicate that the Delayed Detonation scenario can explain 
the major features of observed SNe la, in particular the ex- 
plosion energy and nucleosynthesis (Gamezo et al. 2004}. In 
this scenario, the explosion starts out as a deflagration, with 
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Fig. 12. Equatorial rotational velocity as a function of the mass 
coordinate in the pre-supernova models AAelO, AAe20 and 
AAe30. 



a deflagration-to-detonation transition (DDT) occurring during 
the explosion. However, so far the DDT does not occur selfcon- 
sistently in the explosion models but needs to be put in by hand. 
Lisewski et al. ( 2000 1 showed that a DDT in a SN la may occur 
if turbulent motions with velocities of ~ 10 cm s occur on 
scales of « 10 6 cm. Such high velocities are not expected for 
Rayleigh-Taylor-driven turbulence in SN la deflagration mod- 
els. The rotational velocities in our white dwarf models are 
typically from a few to several times 10 8 cm s _I (Fig. I12> . 
Therefore, even if the white dwarf is spun down by the ex- 
pansion during the explosion, rotation could provide enough 
kinetic energy to enforce the turbulence to the level required 
for a DDT (see also Paper I). 

If the occurrence of a DDT is indeed related to the strength 
of rotation, the DDT might occur earlier (i.e., at higher tran- 
sition densities) in more massive white dwarfs, as those ro- 
tate more rapidly (cf. Fig. I12l> . More massive white dwarfs 
will produce more 56 Ni in such a case, since 1-D models show 
that a higher transition density results in a more energetic ex- 
plosion and thus in the production of more 56 Ni (e.g. Hoflich 
& Khokhlov 1996 1. Even for a pure deflagration, more mas- 
sive white dwarfs are likely to yield more 56 Ni, since they can 
provide more fuel for the nuclear burning. Therefore, we may 
expect the production of more 56 Ni, and thus brighter SNe la 
(Arnett 1982), from more massive white dwarfs. 

The production of more 56 Ni in more massive white dwarfs 
is also required for our super-Chandrasekhar mass scenario to 
satisfy the observational constraints on the kinetic energy (£kin) 
of SN la. Observed ejecta velocities imply E^ n =s 10 51 erg, 
assuming that SNe la come from Chandrasekhar mass white 
dwarfs (e.g. Branch 1992). In other words, SNe la have a spe- 
cific kinetic energy (e^ n ) of about 3.6 x 10 17 erg g" 1 . No clear 
correlation between SN la peak brightness and ejecta velocity 
seems to exist (e.g. Hatano et al. 2000), which might imply that 
SNe la have a similar specific kinetic energy on average. A di- 
versity in the mass of exploding white dwarfs might thus imply 
that more massive white dwarfs yield, in general, more specific 




Fig. 13. (a) Nuclear energy required to explain the specific 
kinetic energy (eid n ) of observed SNe la, as a function of the 
white dwarf mass, i.e, E n = £be + ^wc^kin- The binding en- 
ergy given by Eq. Q2> is used here. Lower and upper bounds 
correspond to ekin = 3.6 x 10 17 erg g _I and 5.4 x 10 17 erg g _1 , 
respectively, (b) Expected nickel mass produced in a SN la ex- 
plosion as a function of the white dwarf mass, assuming that 
Mn; is proportional to E a . Lower and upper limits are com- 
puted from the assumption that minimum and maximum values 
for M Ni at M = 1.38 M are 0.2 M and 0.8 M , respectively. 
The thick solid line in the middle of the gray-colored region 
defines the mean value between these two limits. The dashed 
line gives the mean value of in observed SNe la. Lower 
and upper dotted lines indicate the observationally estimated 
minimum and maximum M^i, respectively (Leibundgut 2000 
see text). 



nuclear energy (e n ) in the explosion, since the specific binding 
energy (cbe) increases with increasing mass (Fig. [9}. Although 
1-D explosion models show no strict correlation between e n 
and the produced mass of 56 Ni (Mm), there is certainly a gen- 
eral tendency that a higher e n corresponds to more M^i (e.g. 
Branch 1992 1. Therefore, larger nickel masses are expected, on 
average, in more massive white dwarfs to explain the observed 
SN la kinetic energies. 

Fig. II 3h shows the required nuclear energy E n to give 
ekin = 3. 6... 5.4 x 10 17 ergg~' as a function of the white 
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dwarf mass, where Eq. fl32i is used for the binding energy 
of white dwarfs at p c = 2 x 10 9 g cm 3 . Here, the adopted 
range of e^ n corresponds to what is typically assumed in the 
literature (e.g. Gamezo et al. 2004 1. The variation of E n in 
this figure remains within 75% in the considered mass range. 
1-D Chandrasekhar mass explosion models by Hoelfich & 
Khokhlov (19961 show that ekin = 3.6... 5.4 x 10 17 erg g" 1 
can be achieved with Mm; ~ 0.2 . . .0.8, within the Delayed 
Detonation scenario. Assuming the average Mm is propor- 
tional to E n , we can estimate a possible range of Mm which 
will yield e kin = 3.6... 5.4 x 10 17 ergg" 1 , as in Fig. [Tib. 
Although the estimated range of Mm is very uncertain and 
needs to be confirmed by detailed explosion models, this con- 
sideration leads us to conclude that rotation has the potential to 
plausibly explain the observed diversity of SN la peak bright- 
nesses: the estimated range of Mm covers most of the observed 
range of 0.1 M Q < Mm < 1.14 M , where the lower and 
upper limits refer to SN 1991bg and SN 1991T, respectively 
(Leibundgut SOOOV 

During the last years, differences in the chemical compo- 
sition - notably, in the C/O ratio or in the overall metallic- 
ity - of exploding white dwarfs have been invoked to explain 
the occurrence of different nickel masses for Chandrasekhar- 
mass white dwarf explosions (Umeda et al. 1999; Hoflich et 
al. |3000; Timmes et al. 2003). However, recent 3-D simula- 
tions by Ropke & Hillebrandt ( 2004 ) indicate that different C/O 
ratios have a negligible effect on Mm- Also the metallicity de- 
pendence appears not able to explain the full observed spread 
(Timmes et al. 125051 . 

Another possible cause for the supernova brightness diver- 
sity within the Chandrasekhar mass scenario might be related 
to differences in the carbon ignition density or in the transition 
density for the DDT. However, for any evolutionary path to- 
ward the Chandrasekhar limit, the accretion rates are restricted 
to 10~ 7 . . . 1O~ 6 M yr _1 , since otherwise a white dwarf can not 
efficiently grow in mass. Since the thermal evolution of non- 
rotation accreting white dwarfs is mainly determined by the 
accretion rate, one can not expect a significant diversity in 
the ignition conditions in Chandrasekhar mass white dwarfs. 
Furthermore, none of these possibilities could explain the large 
amount of M M * 1.0 M in SN 1991T (Fisher et al.[T333}. This 
leaves a diversity of masses above the classical Chandrasekhar 
mass, caused by rotation, as a promising explanation for the 
origin of the diversity in the SN la brightness. 

Also the diversity in the ejecta velocity of SNe la, which 
was attributed to different explosion mechanisms by Hatano 
et al. (2000), can be related to different white dwarf progeni- 
tor masses. Given that E\an - £be _ E n , SNe la with the same 
brightness may have somewhat different amounts of kinetic en- 
ergy, if two different masses result in similar nuclear energy 
and Mni- Similarly, it is possible that some very bright SNe la 
have a smaller kinetic energy than less bright ones, depending 
on how much nuclear energy is released at different masses. 
Therefore, in our scenario, the kinetic energy in SNe la may not 
necessarily correlate with the brightness, in accordance with 
the observations. 

The possible correlation between progenitor mass MsNia 
and the SNe la brightness (Msni 3 - ^SNia correlation) could 



be tested in various ways. Observations indicate that more 
than 70% of SNe la are fairly homogeneous, giving sim- 
ilar brightnesses (Mm — 0.6 M Q ) and light-curves, while 
about 20%. ..30% are SN 199bg-type (sub-luminous) or SN 
1991T-type (over-luminous) (Hillebrandt & Niemeyer. 2000; 
Leibundgut 2000). If the MsNia - ^SNia correlation were to ex- 
plain the observed peak brightness variation, the distribution of 
progenitor masses (Pm SNUi ) should not be flat, but peak at a cer- 
tain mass which is to represent the normal SNe la. Explosions 
from lower and higher masses may then relate to sub-luminous 
and over-luminous SNe la, respectively. Population synthesis 
studies are needed to investigate whether this scenario may re- 
sult in a progenitor mass distribution which can be consistent 
with the observed distribution of nickel masses. 

It is more straightforward to relate the MsNia-^SNia correla- 
tion to the fact that SNe la from early type galaxies are system- 
atically less luminous than those from late type galaxies (e.g. 
Ivanov et al.^000). Since SN la progenitors must be, in gen- 
eral, much older in early type galaxies than in late type ones, 
SN la progenitor systems in early type galaxies should have 
a systematically smaller mass budget, which results in smaller 
A^SNia on average, and thus in less luminous SNe la. Similarly, 
if the progenitor mass distribution changes systematically as a 
function of the metallicity of the progenitor systems, SNe la 
luminosities will also change systematically according to cos- 
mic age. For instance, if the initial mass distribution of white 
dwarfs does not sensitively depend on metallicities, we expect 
the peak of the progenitor mass distribution at a lower mass 
for lower metallicity (Langer et al. 2000). We thus expect dim- 
mer SNe la for lower metallicity, which may have important 
consequences for the cosmological interpretation of supernova 
distances at high redshift. 

The aspherical nature of rapidly rotating white dwarfs 
might, in principle, result in an observable polarization, as 
found in SN 1999by (Howell et al. l200Tl and SN 200 lei (Wang 
et al. I2003i . However, the white dwarf rotation per se may 
not be able to create an aspherical shape of the ejecta: Since 
white dwarf models with p c = 2 x 10 9 g cm~ 3 have a spe- 
cific angular momentum of ~ 10 16 cm 2 s _1 (Fig. |6j, the ro- 
tational velocity of the ejecta will be V e jecta,rot — j/^ejecta — 
0.43 [km s- 1 ](;/3 • 10 16 cm 2 s^/C^ecta/lORo), which is neg- 
ligible compared to the observed ejecta velocities in radial di- 
rection (~ 10 4 kms~'). Instead, an aspherical distribution of 
the nuclear fuel in rotating white dwarfs (see Fig. [2J and a 
consequent aspherical explosion could be more important in 
producing non-radial movements in the ejecta. Future 3-D ex- 
plosion models are needed to relate such effects to our rotating 
pre-supernova models. 

8. Electron-capture induced collapse 

Although the present study focuses on Single Degenerate 
SNe la progenitors, our 2-D models may also facilitate the 
study of double (CO-CO and CO-ONeMg) white dwarf merg- 
ers. During such mergers, dynamically unstable mass transfer 
will result in the formation of a thick accretion disk around the 
more massive white dwarf (e.g. Benz et al. I1990i . While it is 
uncertain how the mass accretion from the thick disk onto the 
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white dwarf proceeds, a rapidly spinning super-Chandrasekhar 
mass merger product can be expected, as recently shown by 
Piersanti et al. (2003b I. In this study, however, the assump- 
tion of rigid rotation again strongly restricts the possible ex- 
plosion masses. Furthermore, the estimated ratio of rotational 
to gravitational energy of ~0.14 in their 1-D rapidly rotat- 
ing white dwarf models exceeds by far the maximum possi- 
ble value for such configurations (Ostriker & Tassoul 1968 , see 
also Sect. 16.4.41 . casting doubts on the general validity of their 
results of the reduced and self-regulated accretion rates due to 
rotation. 

Saio & Nomoto (2004i recently suggested that a double 
CO-white dwarf merger leads to accretion at rates near the 
Eddington limit (~ 10~ 5 M yr _1 ), which, after off-center car- 
bon ignition, will produce a rapidly spinning ONeMg white 
dwarf. When the central density in such an ONeMg white 
dwarf reaches about 4 x 10 9 g crrT 3 , electron capture on 24 Mg 
and 20 Ne begins, and the collapse of the star to form a neu- 
tron star is the expected outcome (Nomoto & Kondo 1991 1. 
Strongly magnetized neutron stars with B ^ 10 15 G and 
P - 1 ms produced by such events are proposed as the cen- 
tral engine of cosmic gamma-ray bursts (e.g. Dar et al. 119921 
Usov lT992l Duncan & Thomson [T992I Yi & Blackman [T998l 
Ruderman et al. 2000 1. While non-rotating ONeMg white 
dwarfs undergo electron-capture (EC) induced collapse when 
their mass reaches the Chandrasekhar-limit, rotating ONeMg 
white dwarfs may need a super-Chandrasekhar mass for such 
collapse. Since a single ONeMg white dwarf may be as massive 
a 1.35M , the formation of ONeMg white dwarfs as massive as 
~ 2.7 M by mergers is, in principle, possible. 

Notably, as shown in Fig. |6] the specific angular momen- 
tum in pre-collapse white dwarfs with Mec.onm ~ 1 -5 M is 
larger than /' a; 1.5 x 10 16 cm 2 s _1 , which is the maximum an- 
gular momentum of a rigidly rotating neutron star. Therefore, 
either a large amount of angular momentum must be lost during 
the collapse, or the resulting neutron stars have to rotate differ- 
entially. In the latter case, the expected mean spin rate ranges 
from a few to several times ~ 10 4 s _1 . 

According to our evolutionary scenarios of accreting CO 
white dwarfs, EC induced collapse is also expected to occur 
if carbon ignition at the center of a super-Chandrasekhar mass 
CO white dwarf is delayed for a long time (Paths B2, C3 and 
D3; see Fig. II 11 1. Fig. [6] shows that the specific angular mo- 
mentum in pre-collapse CO white dwarfs (p c = 10 10 g cm -3 ), 
ranges from 10 16 cm 2 s _1 to 4 x 10 16 cm 2 s _1 where 1.5 < 
M/Mq < 2.0. This is again sufficient for the formation of dif- 
ferentially rotating neutron stars of with spin rates of ~ 10 4 s . 

In both cases, such rapidly rotating, newly formed neu- 
tron stars are believed to lose angular momentum via the r- 
mode instability, magnetic dipole radiation, and/or magnetic 
torques. More importantly, if the differential rotation could 
amplify the magnetic field up to 10 15 G in such newly born 
neutron stars (Ruderman et al. 2000), or sucn fields were 
formed by a dynamo process right after the collapse (Duncan 
& Thomson 1992), the necessary conditions to form a cosmic 
GRBs due to electro-magnetic beaming - i.e., Q as 10 4 s _1 and 
B w 10 15 G- might be satisfied (cf. Wheeler et al. l2000l Zhang 
& Meszaros 2001 1. Although it is a matter of debate whether 



this channel could be a major source for GRBs (e.g. Fryer et 
al. 119991 Middleditch 2004 1, its possibility deserves careful in- 
vestigation. Our 2-D models might provide an ideal starting 
point for such a study. 

9. Conclusion 

We present two-dimensional models of differentially rotating 
massive white dwarfs, with internal rotation laws which are 
derived from evolutionary models of accreting white dwarfs 
(Paper I). In these models, the shear rate inside the white dwarf 
is limited by the threshold value for the onset of the dynam- 
ical shear instability. We find that various physical quantities 
of these 2-D white dwarf models - such as angular momen- 
tum and binding energy - can be accurately expressed as an- 
alytical functions of the white dwarf mass and central den- 
sity, and do not sensitively depend on the details of the rota- 
tion profile (Sect. @}. In particular, we construct models with 
p c = 2 x 10 9 g crrT 3 and 10 10 g crrT 3 , which represent the 
pre-explosion and pre-collapse stages of CO white dwarfs, re- 
spectively, and models with p c = 4 x 10 9 g cm -3 which corre- 
spond to the pre-collapse stage of ONeMg white dwarfs. Based 
on these results, in particular using the mass-angular momen- 
tum relation of our 2-D pre-supernova models, we present new 
scenarios for the evolution of SNe la progenitors (Sect. 16.41 
Fig. II 1> . in the frame of a new paradigm of SN la explosion at 
a variety of super-Chandrasekhar masses. By making use of the 
binding energies of our pre-supernova white dwarf models, we 
suggest that the SN la brightness may be proportional to the ex- 
ploding white dwarf mass on average, and that this correlation 
could explain various aspects of the SNe la diversity (Sect.0. 
We also briefly discussed the electron-capture induced collapse 
of rapidly rotating super-Chandrasekhar mass white dwarf as a 
potential candidate for the central engine of a sub-class of cos- 
mic gamma-ray bursts (Sect. In conclusion, by providing 
2-D white dwarf models at the pre-explosion and pre-collapse 
stages (Sect.[3}, together with analytic relations describing their 
global properties (Sect. |4j, we hope to facilitate future studies 
of such objects, including investigations on the role of the r- 
mode instability and the mass and angular momentum accre- 
tion efficiency for the evolution of accreting white dwarfs, pop- 
ulation synthesis models of the distribution of white dwarf su- 
pernova progenitor masses, its relation to SNe la diversity, and 
most importantly, multi-dimensional simulations of SN la ex- 
plosion and electron-capture induced collapse of rotating white 
dwarfs. 
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Table 1. Properties of selected models in sequence AA. The columns has the following meaning. M: white dwarf mass, J: total 
angular momentum, £2 C : central angular velocity, Q m : moment-of-inertia-weighted mean of angular velocity, W: gravitational 
energy, T: rotational energy, U: internal energy, T/W : ratio of the rotational energy to the gravitational energy, R e : equatorial 
radius, R p /R e : ratio of the polar radius to the equatorial radius. 
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Table 2. Same as in Table^but with sequences ABe, ACe, ADe, AEe and AFe. 
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0.825 


4.042 


1 791 


35.252 


1.567 


96 610 
ZO.OJU 


O OAA 


O AOS 


S16 
U.J JO 


A F7pQ 
ArcV 


1.004 


1.066 


4.620 


zi n 1 7 


37.684 


2.222 


97 ^ 1 Q 

Z / . J Ly 


O O^Q 

u.ujy 


O SA9 
U. J'+Z 


O 478 
U.4 / 


rtrclU 


1 784 


1.431 


5.197 


zL 901 


41.208 


3.201 


98 7zLS 
Zo. / t-O 


O 078 
U.U / 


610 
U.OIU 


O A06 
U.H-UO 


AFp1 1 


1 .yjy 


1.941 


5.775 


A 99A 


45.746 


4.501 


10 9A7 
JU.ZH- / 


O OQ8 
U.U70 


O 68Q 


O 1A8 
U. J^fo 


A Cp 1 1 


9 (\Af. 
Z.U40 


2.318 


6.006 


A 1 QO 

■4-. iyu 


48.912 


5.439 


11 961 
J 1 .ZOJ 


111 
U. 1 1 1 


O 7A1 
U. /Hi 


O 1 1 9 
U. J 1Z 






Sequence AGe: / s h = 




= 0.6 ; fp = 0.01 


; a = 0.9 


, Pc — z* A lu 


gem 






AGel 


1.524 


0.630 


5.775 


3.448 


33.996 


1.121 


96 1Q1 


0.033 


16S 

U.JUJ 


719 

U. / Ji 


Aucj 


1. JOJ 


0.811 


6.237 


3.838 


35.821 


1.606 


97 06Q 
Z / .UOV 


O OzL^ 
U.U4J 


O 1Q6 

u. jyo 


O 6SQ 

u.ojy 


AHp^ 
rWjJCJ 


1 1^7 

i .oj / 


1.032 


6.699 


4.073 


37.862 


2.181 


97 76A 
Z / . / OH- 


O OS8 
U.UJo 


Ai8 

U.H-Jo 


S80 
U. JoU 


a rip a 


i 7^n 


1.269 


6.930 


4.127 


39.786 


2.754 


98 169 
Zo. JOZ 


O 06Q 
U.UO? 


O A86 


O S 1 A 
U.J 1*4 




l.OUJ 


1.527 


7.045 


4.036 


41.535 


3.306 


98 SzLS 
Zo.o'+J 


O OSO 
U.UoU 


O SA1 
U. JH-l 


O AS7 
U.H- J / 




1 SQO 

I.07U 


1.881 


7.161 


3.771 


43.371 


3.934 


9Q 969 
Z7.ZOZ 


O OQ1 


616 
U.OIO 


O 1Q1 
U. jyl 








Sequence Re: Rigid Rotation; p c = 


= 2 X 10 ; 


gem 3 








Kel 


1.3 /O 


0.000 


0.000 


0.000 


29.248 


0.000 


24.JJO 


U.UUU 


U.z9o 


1 nnn 
l.UUU 


Re4 


1.392 


0.180 


1.309 


1.309 


29.715 


0.H8 


24.540 


0.004 


0.306 


0.957 


Re6 


1.415 


0.294 


2.004 


2.004 


30.399 


0.294 


24.802 


0.010 


0.326 


0.884 


Re7 


1.428 


0.344 


2.269 


2.269 


30.765 


0.391 


24.938 


0.013 


0.340 


0.848 


Re8 


1.440 


0.388 


2.472 


2.472 


31.095 


0.480 


25.059 


0.015 


0.357 


0.797 


Re9 


1.449 


0.422 


2.613 


2.613 


31.355 


0.551 


25.151 


0.018 


0.377 


0.754 


RelO 


1.455 


0.443 


2.694 


2.694 


31.519 


0.597 


25.208 


0.019 


0.400 


0.710 


Rell 


1.458 


0.451 


2.724 


2.724 


31.579 


0.615 


25.228 


0.019 


0.426 


0.667 



